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Abstract. - We present a new theoretical analysis of the strongly suppressed F- and M- 
dependent Stark shifts of the Cs ground state hyperfine structure. Our treatment uses third 
order perturbation theory including off-diagonal hyperfine interactions not considered in earlier 
treatments. A numerical evaluation of the perturbation sum using bound states up to n=200 
yields ground state tensor polarizabilities a2(65'i/2, -F) which are in good agreement with ex- 
perimental values, thereby bridging the 40-year-old gap between experiments and theory. We 
have further found that the tensor polarizabilities of the two ground state hyperfine manifolds 
have opposite signs, in disagreement with an earlier derivation. This sign error has a direct 
implication for the precise evaluation of the blackbody radiation shift in primary frequency 
standards. 



Introduction. - Since its discovery, the Stark effect, i.e., the interaction of an atom with 
an external electric field E, has played an important role in the spectroscopic investigation 
of atomic structure. Because of parity conservation, the Stark effect arises only in second 
order perturbation theory for atoms with non-degenerate orbital angular momentum states. 
The energy shift Ai?(7) of a magnetic sublevel I7) = |nLj, F, M) of the hyperfine structure 
is conventionally parametrized in terms of a polarizability 0(7) as 

Ai?(7) = -la(7)E^ (1) 

In second order perturbation theory the polarizability can be decomposed [1] into a scalar 
(2) 

polarizability, , which leads to an F- and M-independent level shift, and a tensor polariz- 
ability, , which yields F- and Af-dependent energy shifts. Because of rotational symmetry, 
the tensor polarizability vanishes for the spherical nSi/2 and nPi/2 states in the alkalis. As 
a consequence, the second order Stark effect in the alkali ground state does not affect its hy- 
perfine splitting, nor does it lift the Zeeman degeneracy of a given hyperfine level. However, 
Haun and Zacharias observed already in 1957 [2] that a static electric field does induce a 
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quadratic Stark shift of the hyperfine transition frequency (-F-dependent Stark effect), and in 
1964, Lipworth and Sandars [3] demonstrated that a static electric field also lifts the Zeeman 
degeneracy within the _F = 4 sublevel manifold of the cesium ground state (an |A/|-dependent 
effect). Improved measurements were performed later by Carrico et al. [4] and Gould et al. [5], 
and were recently confirmed by our own measurements [6-8]. In 1967, Sandars [9] showed that 
the F- and |Af |-dependence of the Stark effect can be explained when the perturbation theory 
is extended to third order after including the hyperfine interaction. His theoretical expressions 
for the third order tensor polarizabilties were evaluated numerically in [3] and [5] under 
simplifying assumptions. This yielded (absolute) values which were systematically larger than 
the experimental values for the five investigated alkali isotopes [5] . It is worth noting that in 
cesium the level shifts due to are almost seven orders of magnitude smaller than the com- 
mon shift of the ground state levels due to the second order scalar polarizability ckp . While 
the scalar Stark shift is understood at a level of (1-2) • 10""^ [10-12], there has so far been 
no satisfactory theoretical description of the tensor polarizabilities, i.e., of the |Af| -dependent 
alterations of the scalar effect. 

In this paper we extend the previous theoretical treatment of the third order effects by 
including both diagonal and off-diagonal hyperfine interactions in the perturbation expansion. 
The numerical evaluation of the perturbation sum then yields a value of the Cs ground state 
tensor polarizability a2 (66*1/2, F — A) which is in very good agreement with the experimental 
results [4-8], thereby bridging the 40-year-old gap [8] between experiments and theory. We 
have also identified an error in Sandars' results [9] concerning the relative sign of the tensor 
polarizabilities of the two ground state hyperfine levels. We discuss the relevance of this sign 
error for the determination of the black-body radiation shift of primary frequency standards 
from measurements of the static Stark shift of the Cs clock frequency. 



Theoretical model. - As shown by Sandars [9], the F- and M-dependent Stark shifts 
can be explained by a third order perturbation treatment, in which the Stark interaction, 
Hgt = — d • E , and the dipole-dipole (HjJ,.) and electric quadrupole (Hjjf) hyperfine interactions, 
are treated on an equal footing. The contributions of Hst and H^f = H^Jf + Hh, can be expressed 
in terms of irreducible tensor operators 

Hst = |e|rC(i).EW (2) 
H^, = (ai^, - yiO [C(2) X S(i)] ] + a,S(i)) • I^^) , (3) 

H^, = 6qC(2).[IWxIW]('\ (4) 

where L, S, and I are vector operators associated with the orbital angular momentum, the 
electronic spin, and the nuclear spin, respectively, and where the C'^'"'^ are normalized spherical 
harmonics of rank k. The first and second terms of Eq. Q represent the magnetic interaction 
of the nuclear magnetic moment with the orbital and electron spin dipole moments respectively 
(dipole-dipole interactions), while the third term represents the Fermi contact interaction and 
Eq. (0J is the electric quadrupole interaction. The corresponding coupling constants are aj.j, 
Os, and hq. The contact interaction term has non-zero matrix elements between S-states 
(L = 0) only, while the first two terms of Eq. Q apply to states with L > 0. The electric 
quadrupole term of Eq. Q requires L > 1 and J > 1, and for Cs states its matrix 

elements are two orders of magnitude smaller than all other contributions. 
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The third order energy perturbation of the sublevel \a) — |6S'i/2, M) is given by 

where Ep and i?-^ are the unperturbed state energies. Of all the terms obtained by substituting 
W = Hgt + Hi,f into Eq. only those proportional to give nonzero contributions to the 
Stark interaction, and the sums have to be carried out over all states according to the selection 
rules imposed by the operators. 

We first address the second term of Eq. (jsj , whose only non- vanishing contribution is 

where |/9) = \nPj, f, m). Diagram A of Fig. ^ shows a graphical representation of this term in 
which the Fermi ground-state contact interaction appears as a multiplicative factor, making 
this contribution i^-dependent. Comparing the sum in Eq. © to the expression 

^^(2)(^)^^MHs^, (7) 

for the second order scalar Stark effect, one sees that the F-dependent third order shift is 
suppressed by a factor on the order of Efif{6S)/{E(iP — E^s) ~ 10^^. Following the definition 
of Eq. 1^ the contribution of Eq. ((HJ can be parametrized in terms of an i^-dependent third 
order scalar polarizability (651/2, E). This F-dependent term gives the main contribution 
to the Stark shift of the hyperfine transition frequency. 

We next address the first term of Eq. jSJ- The nonvanishing contributions involving 
diagonal matrix elements of H,,f are given by 

2^ {P\ Hm IP) -j^ TTT^ = 2^ Ehf{f3)— — ^ , (8) 

where = \nPj, f,m). They are represented by diagram B of Fig. ^ and are suppressed by 
a factor on the order of Ehf{6P)/{EQP — Eqs) ~ 10~^ with respect to the second order scalar 
shifts. The electric field dependent factor of Eq. ((Sj) has only a rank fc = (scalar) contribution, 
while the dipole-dipole [Eq. (0)] and the electric quadrupole [Eq. parts of the factor 
in Eq. JH)) have the rotational symmetries of A; = 0, 2 and k = 2 tensors, respectively. The 
contraction of the hyperfine and the Stark interactions in Eq. ijHJl thus yields both scalar and 
second rank tensor contributions. The scalar part of Eq. ||SJ) has the same -F-dependence [8] as 
aQ {E), and gives a correction to the latter on the order of 1%, while the second rank tensor 
part has an F- and A/-dependence, which can be parametrized in terms of a third order tensor 
polarizability aj'^^ (65*1/2, F). 

The total third order polarizability can thus be written as 



.(3) (65i/2, F, M) = (65i/2, F) + ai'^ (65i/2, P f^lj^f/^^^^^ 



where the dependence on the angle 9 between the electric field and the quantization axis is 
given by f{0) = 3cos^ 9—1. Equations H6I8|I can be reduced by applying the Wigner-Eckart 
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theorem and angular momentum decoupling rules. For cesium (/ = 7/2) the explicit F- and 
M-dependences of the third order polarizabilities for 9 = are 

'^M^ — 90 

a^^\F = A,M) = ao + (ai+a2) , (10a) 

J3),p_. 9 / 5 \3M2-12 



a^'\F = i,M) = --ao+ (^-ai + -a2j . (10b) 

The scalar part can be expressed in terms of radial integrals 
■J ^ 

"0 = ^ X! [^^"^1/2 (3^6Si/2 + ^"^1/2) + 2C„P3/, (3^651/2 - 5A„P3^ J] , (11) 

n 

where ^6Si/2 ^-nd Anp. are the hyperfine coupling constants and where 

_ |^6S,nPj I 



(12) 



■Res.nPj is the radial integral between the ground state and the excited \nPj) state. The Fermi- 
contact interaction (proportional to ^651/2) provides the dominant contribution to oq which, 
as mentioned, also has a small contribution (proportional to Anp^ ) from the scalar part of the 
magnetic dipole interactions. Eqs. © and ijSJ) can be expressed in a similar way as 

ai = -— ^ [(2^nPi/2CnPi/2 - SAiPj/aCnPa/s) + {2AnP^^^CnP^^^ + AnP^^^CnP:,/^)] , (13) 
n 

where we have explicitely separated the contribution of the orbital part of the interaction 
(first term) from that of the spin dipolar part (second term), and 



02 



- Q X! ^nP3/2<^nP3/2 ' (14) 



where Bnp^^^ ^^'^ quadrupole hyperfine coupling constant. 

Equations H10() closely resemble the expressions derived by Sandars [9], except for the 
negative sign of the oi term in Eq. (jlOb(l which is positive in Sandars' work. To our knowl- 
edge, this sign cannot be derived from any prior experiment. We have recently confirmed 
experimentally [7] that the tensor polarizabilities of the F = 3 and F = 4 states have indeed 
opposite signs as predicted by our calculation. This sign error has remained unnoticed in the 
literature for almost 40 years and we will come back to its relevance for atomic clocks below. 

Numerical evaluation of the tensor polarizability. - The tensor polarizability (oi term) 
of the F = 4 state of cesium was evaluated in [5] by considering only diagonal matrix elements 
of for the 65*1/2 and the 6Pj states. The contribution of the orbital magnetic dipole 
hyperfine interaction [first term in Eq. H13|) ] was neglected, as well as the spin-orbit splitting 
in the denominators of Eq. ((Sj. The authors assumed furthermore that Aqp^^^/Aqp^^^ = 5, 
valid for one-electron atoms, while for Cs the corresponding ratio of experimental values is 5.8. 

(2) 

Under the latter two approximations one can factor the second order scalar polarizability 
out of the expression for the tensor polarizability. Those simplifying assumptions yielded the 
a'^^ (F = 4) value represented as point (f ) in Fig. |21 in disagreement with the experimental 
results. We have reevaluated the result by including the orbital part of the hyperfine inter- 
action and then rescaling the value of [5] by using the most recent experimental value of 
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Fig. 1 - Contributions of diagonal (A and B) and ofT-diagonal (1-5) hyperfine matrix elements to the 
third order Stark effect. Dotted/solid lines represent the Stark and hyperfine interactions respectively. 
The relative contributions of the diagrams are given in the text. 



[10]. This yields the value (f) in Fig. [21 thus increasing the gap between theory and ex- 
periments. In a second, more precise, calculation we dropped all the simplifying assumptions 
mentioned above [8], still keeping diagonal matrix elements only, and numerically evaluated 
Eqs. H13I14II by using recent experimental values [13] of the reduced dipole matrix elements 
(65*1/2 II f^ll As a result, the discrepancy with experiments becomes even larger [point 

(f") in Fig. 12], and does not change significantly when extending the perturbation sum to nPj 
states with n > 6. 

Off- diagonal hyperfine interaction. - All calculations described above considered only 
diagonal hyperfine matrix elements. However, the first term in Eq. (O allows also off-diagonal 
hyperfine terms. FigureQlgives a schematic overview of all possible off-diagonal configurations 
compatible with the hyperfine and Stark operator selection rules. It is interesting to note that 
some off-diagonal hyperfine matrix elements (diagrams 1 and 2, Fig.^ were already considered 
by Feichtner et al. [14] in their calculation of the clock transition Stark shift, but for unknown 
reasons such terms were never considered in the tensor polarizability calculation. 

The contact interaction selection rule Ai — admits only off-diagonal matrix elements 
between the ground state and higher lying 5*1/2 states (diagram 1, Fig.^l. The orbital and spin 
dipolar magnetic interactions also obey AL — and have nonvanishing off-diagonal matrix 
elements of the form (/3i| H^f |/32) where 1/3^) = \niPj.) (diagrams 2, 3 and 4). Due to the second 
rank tensor character of C^^^ the spin dipolar term also couples states with AL = ±2 and 
can thus contribute to the third order Stark effect with off-diagonal matrix elements between 
the ground state and I?3/2 states (diagram 5). We have used the Schrodinger equation with a 
statistical Thomas- Fermi model potential to calculate the relevant wave functions of the free 
Cs atom. Corrections for the dipolar and quadrupolar core polarization as well as spin-orbit 
interaction with a relativistic correction factor following [15] were included. The electric dipole 
and hyperfine matrix elements were calculated using the Schrodinger wave functions for all 
matrix elements for which no precise values could be found in the literature. In this way we 
evaluated all diagrams in Fig. ^ by running the summation indices m and n up to 200. 
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Relative contributions and numerical results. - Diagrams A and 1 of Fig. ^ yield only 
F-dependent energy shifts, and thus do not contribute to the tensor polarizability. All other 
diagrams produce F- and ill-dependent effects. The relative importance with which diagrams 
B and 2-5 contribute to a':^^ is +145, +99, -40, -3, and -f Of %, respectively. A numerical eval- 
uation of the perturbation sum with all diagonal and off-diagonal matrix elements mentioned 
above gives 

= 4) -3.72(25) x fO"^ Hz/(kV/cm)2 , (15) 

for the tensor polarizability, in which the contribution of 02 is 2 x 10^^ Hz/(kV/cm)^. This 
result is shown as a dotted line in Fig. [5] together with previous theoretical and experimental 
results. We estimate the uncertainty of our calculated value to be 7%, based on the precision 
(2-8%) with which our Schrodinger solutions can reproduce measured dipole matrix elements 
and hyperfine splittings, and considering that some (more precise) experimental values were 
included in the calculations, whenever they were available. We have also verified by explicit 
calculations using continuum wave functions that continuum states contribute only at a level 
of lO^"* to the diagrams relevant for aj^-* (~ ai), while it was recently shown [16, 17] that the 
continuum contributes significantly (« 10%) to diagram 1, which dominates the BBR shift 



(3) 



ao) term in Eq. ^ [18]. 



via the a, 

The corrected sign of the ai terms has an important implication for frequency standards. 
From Eqs. (|10|l the static Stark shift of the Cs clock frequency is given by 



Auoo = Au (65i/2, 4, 0) - Au (65i/2, 3, 0) 



16 4 ■ 
y ao - y ai / [0) 



(16) 



One of the leading systematic shifts of the clock transition frequency is due to the interaction of 
the atoms with the dynamic Stark interaction induced by the blackbody radiation (BBR) field. 
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Fig. 2 - The third order tensor polarizability of the _F = 4 Cs ground state. The filled squares 
(■) represent experimental values of Carrico et al. [4] (a), Gould et al. [5](b), Ospelkaus et al. [6](c) 
and Ulzega et al. [7](d). The empty square (□) (e) represents their weighted average. The dots (•) 
represent the theoretical value from [5](f), and our reevaluations (f, f") discussed in the text. The 
dotted horizontal line is the result of the present work with its uncertainty (shaded band). 
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This shift can be calculated from Eq. (|16|) by averaging E'^ over the Planck spectrum [16,17, 19] . 
Because of the isotropy of the blackbody radiation the dependence in Eq. \i<6\ vanishes, so 
that the BBR shift is determined by the scalar part only. The BBR shift coefficient can be 
deduced from the measured shift of the clock transition frequency induced by a static electric 
field. For this one has to know the value of the ai term in Eq. which is affected by the 
sign error, the coefficient of ai being — ^ when derived from Sandars' formula and — | from 
our calculation. The most precise measurement of the static Stark shift was performed by 
Simon et al. [20] from which the authors extracted — |ao — —2.273(4) kHz/{kV/cm)'^, while 
our evaluation with the correct sign yields — |ao = —2.281(4) kHz/{kV/cm)'^. The correction 
of the sign error thus changes the BBR shift rate by a value which is twice as large as the 
reported experimental uncertainty. 

Conclusions. - We have extended a previous treatment of the Stark interaction by 
including off-diagonal hyperfine matrix elements in the third order perturbation expansion. 
Our calculation of the tensor polarizability yields a good agreement with all experimental 
data, thereby removing a 40-year-old discrepancy. A sign error identified in Sandars' model 
leads to a new expression for the static Stark shift of the hyperfine transition frequency which 
requires Cs primary frequency standards to be corrected at a level of 6 x 10"^^, which is below 
their present accuracy. 
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